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I. INTRODUCTION 



The constraints imposed by chiral symmetry breaking in gauge theories can be surprisingly strong. Low-energy 
theorems, the dynamics of pseudo-Goldstone bosons in an expansion around the zero-momentum limit, and the whole 
framework of effective chiral Lagrangians are examples of this. Generally these constraints are imposed on the effective 
low-energy degrees of freedom only. It is much more surprising that both spontaneous chiral symmetry breaking and 
the explicit breaking of chiral symmetry due to the U(l) anomaly also can be used to give exact analytical predictions 
for the underlying fermion degrees of freedom. This is possible when one restricts the gauge theory to a large but 
finite four-volume V obeying the inequality |]J 

V « -\ , (1) 

where is the pseudo-Goldstone mass. In this rather extreme limit the QCD partition function depends on the 
fermion masses m t only in the particular combination fii = VEuii, where S is the chiral condensate. While the 
four-volume V must be taken to infinity in order to obtain analytical predictions, a finite-size scaling regime is thus 
achieved by sending fermion masses rrii to zero at just such a rate that the /i^'s remain fixed. This is an exact finite- 
size scaling region in the same sense as near critical points: we can reach as accurate agreement we wish by simply 
choosing a sufficiently large four-volume V. In contrast to what one is accustomed to in statistical mechanics, the 
finite-size scaling functions are in this case known exactly, both in shape and absolute normalization, once one has 
the value of the infinite- volume chiral condensate S. 

One of the remarkable aspects of the finite-size scaling region (1) is its relation to Random Matrix Theory results 
that have proven to be universal |§-0|. There is a beautiful classification of the possible chiral symmetry breaking 
patterns for different gauge groups and color representations of the fermions in terms of the classical Random Matrix 
Theory ensembles labeled by the so-called Dyson index. This leads to three major universality classes ||. For 
the quenched case, the analytical prediction for the mass-dependent chiral condensate was in fact first obtained by 
Verbaarschot H using the exact formula for the microscopic Dirac operator spectrum as derived from Random Matrix 
Theory. It has later become clear that these results can also be derived from finite-volume partition functions alone 
. A powerful analytical technique uses fully or partially quenched (supersymmetric) chiral Lagrangians Q . Lattice 
gauge simulations have already shown nice agreement with the exact analytical predictions for the microscopic Dirac 
operator spectrum associated with all three different universality classes f^-flll using staggered fermions. It has 
been particularly challenging to see also the detailed analytical predictions in gauge field sectors of fixed non-zero 
topological charge v, and this has very recently been achieved using perfect actions and overlap fermions 

The purpose of this paper is to perform a systematic series of lattice gauge theory tests of the exact predictions 
related to the mass-dependent chiral condensates and one of the chiral susceptibilities. For the case of gauge group 
SU(3) and staggered fermions in the fundamental representation such an analysis was first performed || on the basis of 
lattice gauge theory data from the Columbia group . These data, based on configurations with Nf = 2 dynamical 
fermions, were taken for finite-temperature lattice volumes, and with rather large correlation lengths. Because the 
dynamical fermion masses were large on that scale, and the physical pions therefore not obeying the inequality (1), the 
configurations were in fact to be considered as completely quenched on the scale of the mass of the valence fermions. 
As there are very accurate estimates for the infinite-volume chiral condensate £ for the same theory at different j3- 
values JTof , we now have parameter-free predictions at these (i- values. Furthermore, we can probe much larger physical 
volumes (and hence obtain higher accuracy), and our result will not be contaminated by finite-temperature effects. 
This case corresponds to the chiral unitary (chUE) Random Matrix Theory ensemble. We next turn to gauge group 
SU(2) and staggered fermions in the fundamental and adjoint representations. The former case corresponds at our 
finite lattice spacings to the chiral symplectic ensemble (chSE) in the Random Matrix Theory classification, while the 
latter corresponds to the chiral orthogonal ensemble (chOE). In this way we cover all three major universality classes. 
However, staggered fermions suffer from two significant defects in this context. First, at our finite lattice spacings the 
staggered fermions with SU(2) gauge group do not fall into the right universality classes as compared with fermions 
in the continuum ||. Second, artifacts due to finite lattice spacings prevent us from testing more than the gauge field 
sector of topological charge v = with these fermions. Both of these shortcomings can be overcome by the use of more 
sophisticated fermion formulations. We shall here provide lattice data obtained with overlap fermions jl5| . Here finite 
lattice-spacing analogs of continuum relations for the chiral condensate in non-trivial topological backgrounds can be 
established [ |l6| . We thus simultaneously achieve both the correct identification with respect to continuum universality 
classes, and correct relationships in non-trivial topological gauge field sectors. We shall throughout restrict ourselves 
to the quenched limit, Nf — 0. Analytically this limit is readily taken, both from Random Matrix Theory and from 
the (quenched) finite-volume partition functions, and the answers have been shown to agree. There are thus precise 
and unequivocal analytical predictions also for this case. 
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For which observables do we have exact analytical predictions? Essentially all quantities for which the partition 
function itself Z({pi}), perhaps extended with additional fermion species, is a generating function. The simplest 
quantity to focus on is obviously the mass-dependent chiral condensate itself, which in the quenched theory simply 
reads 

%U|m[Z(,)]. (2) 

Here £ is the genuine chiral condensate in the massless limit of the infinite-volume theory, and p = VT,m is the 
microscopically rescaled quenched "valence" fermion mass m. Because of the relation (at fixed topological charge v) 



W d C ^ + ^, (3) 



the mass-dependent chiral condensate tests a massive spectral sum rule for the microscopic density p^ (C; p) of 
the Dirac operator spectrum [^],[l7]] .[] We shall here restrict ourselves to the quenched cases, where pi" is mass- 
independent (and the partially quenched cases are completely analogous). This is one way in which the analytical 
prediction for the quenched limit can be obtained (the other proceeds directly from the quenched chiral Lagrangian 
P|) . In numerical simulations the condensate measurements are of course performed very easily (from the trace of the 
propagator), without ever having to compute the Dirac operator spectrum itself. Nevertheless, we have sometimes 
found it convenient to supplement direct measurements by appropriate sums over eigenvalues. In this connection it 
is important to stress the following point. Because of Eq. (|), we are in fact concerned with tests of the analytical 
predictions for p s (£,ii). The advantage of using S„(/i) to test these predictions is that it in a quantitative manner 
probes the microscopic spectral density in different regimes. For instance, by going to very small values of p the 
condensate becomes very sensitive to the precise eigenvalue distribution around £ ~ p. In particular for the chUE (in 
the v = sector) and the chOE (in the v = and v = 1 sectors) £„(/z), as we shall see below, is extremely sensitive 
to the low-£ distribution of the smallest non-zero eigenvalues. This effect can be enhanced by considering in addition 
an observable such as a chiral susceptibility, as we shall discuss below. 

Considering the expression (|^) for the chiral condensate, one might wonder about the necessity of subtractions. 
After all, even in free field theory the spectral density of the Dirac operator goes like p(X) ~ A 3 , and the spectral 
representation of the condensate is thus ultraviolet divergent. There are no such divergences in the finite- volume 
scaling regime considered here, and we should make no subtractions in the chiral condensate either. Although this 
point was already explained in the paper of Leutwyler and Smilga [jlj, it is worthwhile to repeat it here. The 
explanation is as follows: What we are computing here is not the conventionally defined chiral condensate. We are 
taking a correlated limit of V — > oo and m — > such that p = mY<V is kept fixed. In this limit the condensate £(/i) 
receives contributions only from Dirac operator eigenvalues on the scale of m <C Aqcd, and below. The ultraviolet 
end of the Dirac operator spectrum is not ignored: the corrections (and hence subtractions) from this region are of 
the kind mA 2 and m 3 In A, where A is the ultraviolet cutoff. In the scaling region where p — m£ V is kept fixed, 
these terms are suppressed by 1/V and 1/V 3 , respectively. In other words, the ultraviolet end of the Dirac operator 
spectrum in this region enters only as 1/V corrections to the main predictions. In addition, there are of course also 
1/V corrections from the neglect of non-static modes in the effective partition function, so all of these 1/V corrections 
are effectively beyond our control. Consistent with this observation is the fact that 1/V corrections are non-universal 
in the Random Matrix Theory context ]l8[| . Here we are interested only in the leading, universal, predictions for 
V^oo. 

Corresponding to the three different universality classes, there are three distinct predictions for the mass-dependent 
condensate. We shall here give the predictions for gauge field sectors of fixed topological index v, and for any number 
of massless flavors Nr. First, for gauge groups SU(iV c ) with N c > 3 and fermions in the fundamental representation 
the prediction reads || 

v = M [lN f +u{p)K Nf+v (p) + I Nf+ v + i{p)K Nf+v -i 

(/*)] + 7 . ( 4 ) 

Zj p 

where I n (x) and K n (x ) are the two modified Bessel functions. This is the universality class of the chUE in the Random 
Matrix Theory classification. Staggered fermions in the same representation and for the same gauge groups are here 



1 In what follows we shall for convenience, to avoid absolute-value signs, always consider v non-negative. 
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belonging to the same universality class as continuum fermions. In this case the explicit form of the microscopic 
spectral density pa {CifH) °f a l so Nf massive fermions (of masses pi) is known, so that one has also available 
complete analytical predictions for partially quenched chiral condensates with massive fermions. These can also be 
derived directly from partially quenched chiral Lagrangians (see ref. @]). 

For the chSE universality class, where the microscopic spectral density for Nf massless fermions in a sector of 
arbitrary topological charge v has been given in very compact form in U, we have been able to reduce the chiral 
condensate to the following. When Nf is even we find 



^ChSE 



(m) 



= 2p [l Nf+2v (2p)K Nf+ 2v(2p) + I Nf+2u+ i(2fi)K Nf+ 2 V -i(2fi)] + - 

P 

- 2(-iy +N f/ 2 K Nf+2 ,A2p) [p (l {2p) - |(L (2 M )/i(2//) - Li(2 M )I (2 M )) 
- £ (-l) k hk+i(2p) 

k=0 



where L„(a;) denotes the nth order modified Struve function. For Nf odd we find 

£f SE ( u ) 



(5) 



2/1 [l Nf+2 i,(2fl)K Nf+ 2i,(2fJ,) + I Nf+ 2v+l(2p)KN f +2v-l(2p)] + - 

u+(N f -X)/2 

l-I (2p)-2 i-^) k hk(2p) 



(-ir + ^ +i ^ 2 K Nf+2u (2p) 



(6) 



These predictions pertain to gauge group SU(iV c ) with N c > 2 and fermions in the adjoint representation (for contin- 
uum fermions). This is also the universality class relevant for staggered fermions in the fundamental representation 
and gauge group SU(2). Finally, the universality class of the chOE predicts a chiral condensate of the following form. 
For v odd we find, using the recent compact expression for the microscopic spectral density of that case [jl9| , 



^ChOE 



(P) 



= P [l2N f +u(p)K 2 N f +u{p) + l2Nf+L>+l(p)K2Nf+v-l(p)] + 



+ (-1)"' 



P 



N t +(v-l)/2 

I 2 N } +Ap)Ko{p) + 2 {-^) k hN s +Ap)K2k{p) 

k=l 



(7) 



while for v even we have been able to reduce the answer to 
^ OE (p) 



= P [l2N f +u(p)K 2 N f +u{p) + l2Nf+f+l(p)K2Nf+v-l(p)] + 



N f +v/2-l 



,2fc+l 



fe=0 



(2N f 

{-l) Nf+lj/2 K 2Nj+v {p) [p (l {p) - ^(L (p)h(p) - Li^/oM) 

Nf+v/2-l 

-2 (-V k h k+ i{p) 

k=0 



(8) 



where, with the usual convention, (—1)!! = 1. Eqs. (Q) and (||) give the chiral condensate in SU(2) gauge theory and 
continuum fermions in the fundamental (pseudo-real) representation. They also correspond to gauge group SU(iV c ) 
with N c > 2 and staggered fermions in the adjoint representation. 
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FIG. 1. The condensate for staggered fermions in (a) the fundamental representation of SU(3) as function of the fermion 
mass, (b) rescaled form as a function of (i = m'EV, (c) the fundamental representation of SU(2), (d) the adjoint representation 
ofSU(2). 



4 



II. STAGGERED FERMIONS 



Analytical predictions for the quenched chiral condensate and higher chiral susceptibilities are all restricted to 
sectors of fixed gauge field topology.0 As mentioned above, with staggered fermions lattice simulations of the chiral 
condensate at realistic values of the coupling do not see any trace of gauge field sectors except the one of topological 
charge v = [0]. The comparisons one can make are therefore slightly limited. On the other hand, computationally 
the staggered fermion formulation is extremely convenient for our purposes. We shall therefore start with a systematic 
study of the chiral condensate based on staggered fermions. 

For the determination of the staggered quenched chiral condensate, we used a stochastic estimate method for several 
3 in SU(2) and SU(3) with fermions in the fundamental and adjoint representations. We are interested in solving 
the linear system of equations (m)D(m)ri(m) = (Z)t(0)Z)(0) + m 2 ) 77(771) — b for some stochastic source b and the 
staggered Dirac operator D. We can relate the solution 77(771) in a fixed background gauge field to the final quantities 
we are interested in with the following expressions: 

(&|Z? _1 (m)|6) = m b^r/irn) 
{b\(D^m)D{m))- 1 + D- 2 {m)\b) = 2 m 2 7^(777)77(777) (9) 

A multi-shift conjugate method |20) was used to solve the required linear system for several fermion masses. We 
supplemented these measurements with computations of the lowest 50 eigenvalues of the staggered Dirac operator, 
typically in 32 bit single precision. These eigenvalues were used to compute a truncated spectral sum approximation 
for the condensates at small fermion masses. The truncated spectral sum method greatly improved the statistics. 
However, for the small lattices (4 4 and 6 4 ) typically on the order of 100000 configurations were needed. As will be 
discussed below, this reasonably large amount of statistics was necessary to adequately sample the lowest eigenvalue 
distribution of the Dirac operator probed by the small fermion masses used for our tests of the predictions of Chiral 
Random Matrix Theory. In addition, for the orthogonal ensemble cases, double precision was used for the eigenvalues. 

Consider first the gauge group SU(3) and staggered fermions in the fundamental representation. The universality 
class is thus that of the chUE, the same as in the continuum. In Fig. we show the quenched chiral condensate 
as a function of the valence mass m for a few different /3-values, and various different lattice sizes. At the shown 
values of the lattice coupling 8 the infinite- volume chiral condensate S has already been determined to high accuracy 
from independent studies [jlCj] . This means that the finite-size scaling function Sj,(/i) of Eq. (^) is parameter- free. 
The first observation is that all the different data roughly collapse down on one universal scaling curve, once plotted 
against /1 = mTiV as in Fig. [ljj. And the analytical prediction (Q) for this curve is remarkably well reproduced by our 
lattice data. This of course just confirms on these lattice volumes and these /3-values the observation first made by 
Verbaarschot || on the basis of data from the Columbia group. It is particularly interesting to look at the small-mass 
behavior. If we expand the condensate in (|J) for small /i, we get (for v = 0) 

^^ = -(M§)+7-i)M + 0(M 3 ), (10) 

where 7 is Euler's constant. There is the expected term linear in 77, but in addition a logarithmic correction of the 
form 7_iln(/i). This latter term, which should not be confused with so-called "quenched chiral logarithms", arises from 
the infrared part of the integral in Eq. (||), and is thus very sensitive to the fall-off of p s (C) near £ ~ 0. Indeed, for 
the chUE the quenched microscopic spectral density reads Q : 

p1 0) (C) = ^[Jo(C) 2 + Ji(0 2 ] , (11) 

which for small values of C, behaves like p^(C) = C/2 — C 3 /8 + - ■■• The leading linear term here is responsible for the 
/7ln(/_i) piece in ([To]). The chUE cases with v > 1 lend themselves more easily to an analysis of the low- mass behavior 
of the condensate, since the integrals 



2 With dynamical fermions one can analytically perform the required sum over topology 111 , but this is not possible in the 
quenched case without additional assumptions about the distribution of winding numbers [|L8| . 



5 



0.6 




0.4 



c 

1 

Oh 



0.2 



0.0 



0.00 0.01 0.02 0.03 

ASV 



0.04 



0.05 



FIG. 2. The distribution of the lowest eigenvalue for staggered fermions in the adjoint representation of SU(2) at ft = 2.0 
and 4 4 corresponding to the orthogonal ensemble. The scale is greatly enlarged to show the leading edge of the distribution. 



in those cases are convergent. We can thus make the following rewriting: 
S Jo C + A 1 A* Jo 



1 ^ 



2 



C 2 C 2 (C 2 +m 2 ) 



- , (12) 
A* 



and consider each term separately. It follows that the leading, linear, term in the expansion of E„(/i) for v > 1 has as 
coefficient the first Leutwyler-Smilga sum rule (as extended to this quenched case) |l]||: 



^ hUE (A0 




(13) 



We next turn to gauge group SU(2) and staggered fermions in the fundamental representation. The analytical 
prediction is here that of the chSE universality class, with a chiral condensate as in Eqs. (§) and (§). Wc again choose 
/3-values for which the infinite- volume chiral condensate £ is known to high accuracy, so that the analytical predictions 
(^) and (^) also are parameter-free. One remark is in order here: in the chSE universality class every eigenvalue is 
doubly degenerate. The analytical predictions from RMT consider only one of the eigenvalues from each degenerate 
pair. A stochastic estimate of the condensate, on the other hand, contains the contribution from both eigenvalues of 
each pair, and is thus a factor two larger. We have therefore divided the stochastic estimate of the condensate by this 
factor two in order to compare to the analytical prediction. In Fig. ^c we show how all data nicely collapse down on 
the universal scaling function (|^) for Nf = 0. The agreement is seen to be extraordinarily good over more than three 
orders of magnitude. 

The microscopic spectral density of that case Q reads as follows for Nf = 0: 

pf\0 = CJi(2C) 2 - y (H (2C)Jo(2C)Ji(2C) - H 1 (2C)J (2C) 2 ) , (14) 

where H n (x) is the nth order Struve function. Since the small-C expansion is of the form p^\C) = C 3 /3 + • • ■, we 
can make the same rewriting as above (see Eq. ([l2])) and consider each term separately. It follows that also here the 
leading, linear, term in the expansion of So(/x) has as coefficient the first Leutwyler-Smilga sum rule for that case 
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FIG. 3. The quenched susceptibility ui for staggered fermions in (a) the fundamental representation of SU(3) as function of 
the fermion mass, (b) in rescaled form as a function of \i = mT.V , (c) the fundamental representation of SU(2), (d) the adjoint 
representation of SU(2). 
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(15) 



This linear behavior with coefficient one is precisely what is observed in Fig. ^p. The same argument goes through 
also for sectors of non-trivial winding numbers v , in which case the formula reads: 




A* l + 2M r 

We finally present lattice gauge theory data for the SU(2) gauge group, and staggered fermions in the adjoint 
representation. Here data should align on the universal scaling curve of the chOE universality class (see Eq. (H)), and 
we show the results of a few high statistics (but rather small) lattice volumes and two different (3- values in Fig. [Lp. The 
/3-values were again chosen on the basis of having already a good estimate for the infinite-volume chiral condensate 
S (for the adjoint representation) H]. The analytical curve is seen to have a surprising behavior: it rises, even here 
in the v = case, with decreasing (rescaled) fermion mass \i. This unusual feature is a reflection of a peculiarity of 

the quenched microscopic spectral density p a (C) f° r the chOE (see the third of ref. p| and e.g. p9|), 



^H0 = ^M0 2 + \m0 



1-^(^(0^(0-^(0^(0) 



(17) 



Contrary to all other microscopic spectral densities for the chiral ensembles, the above function does not vanish at 
C = 0. This is an artifact of the quenched limit, and it implies that in this quenched theory one can have spontaneous 
symmetry breaking even if one is taking the limit V — > oo and m — > in a correlated manner. Conventionally 
the possibility of spontaneous symmetry breaking implies that one first sends the volume V to infinity, and only 
subsequently takes the massless limit, m — > 0: 

S = lim lim (tpip) . (18) 

rn — >0 V — >oo 

In the quenched case corresponding to the chOE we observe that spontaneous symmetry breaking can occur even if 
we take the simultaneous limit V — > oo and m — > 0, with mV fixed. This holds only in the v — sector. For v ^ 
we face the usual situation that the chiral condensate diverges like I/a*. This holds in the quenched theory as well 
when we sum over topological charges Q. 

Also in this case we can analyze the limit of \x — > analytically for v = 0. The reason for the unusual phenomenon 
of a constant mass-dependent chiral condensate in the limit — ^ is the term Jo (0/2 in Eq. (|l7]). It is this term 
that leads to a non-vanishing microscopic spectral density at Q = 0, and one can easily confirm that it is also this 
term that is responsible for the leading small-At behavior of the chiral condensate in this case. Using 



^^^2 = ^(Io(/M-L„(/0) (I')! 



and the small- a* expansion of the modified Struve function L (/i) = 2A«/7r + . . ., we see that only the first piece 
contributes in the limit fj, — > 0. From Eqs. (g) and ([l7|) we finally get^. 

= £ + (20) 



for this universality class. An approach towards this constant value is seen in the data of Fig. ^d, but the signal 
obviously gets rather noisy around fi ~ 10~ 3 for these lattice volumes. We show in Fig. || magnified plot of the 
distribution of the smallest eigenvalue along with the curve for the fit for the 4 4 lattice. We see there is a reasonable 
sampling of the distribution for the very small eigenvalues, but even larger statistics beyond our 135000 configurations 
are needed to really adequately sample this region and hence give very reliable estimates for the condensate. Again 
we see that the chiral condensate is an extremely sensitive probe of the smallest Dirac eigenvalue spectrum. For 



3 Note that the fi — > limit gives a condensate that is a factor of tt/2 larger than the conventionally defined chiral condensate. If 
finite- volume effects do not eventually cut off the lowest eigenvalue in this case, one can even have spontaneous chiral symmetry 
breaking without first taking the infinite- volume limit (we thank J. Verbaarschot for emphasizing this last point). 
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example, the statistical fluctuation that causes a small surplus of eigenvalues very close to the origin in Fig. || reflects 
itself directly in the slightly larger chiral condensate in Fig. |l|d. The deviation is seen clearly on the non-logarithmic 
vertical scale. 

In sectors of non-vanishing topological charge v, the microscopic spectral density vanishes at the origin, and if 
it were not for the ^//i-piece, the mass-dependent chiral condensate would then also vanish as /i — ► 0, even in the 
infinite-volume limit. For example, for v = 1 the expansion for small [i reads in this case 



1 



m(f)+7 



(21) 



with, again, a /i ln(/i)-term in addition to the purely linear contribution. 

Other physical observables can of course be extracted from the finite- volume partition function. We shall here focus 
on one such observable, a chiral susceptibility w(/i) which we define as: 



Ps 



(C;aO 



(C 2 + m 2 ) 2 



2v 
7 



(22) 



We expect this quantity to be a more sensitive probe of the rescaled eigenvalues £ at a specific rescaled mass /j, than 
the chiral condensate because of the higher power occurring in the denominator of the integrand. This quantity is 
especially easy to compute in the quenched limit, where the spectral density is /i-independent. One then has 













d 




^ dp. 




"X^Gu) 




A 4 



d(- 



o (C 2 +M 2 ) 

- KM 



(23) 



That combination is particularly useful in testing the small deviation from linear behavior of X y (/i) in, for instance, 
the case corresponding to the chUE with v = 0. In general, for a quenched condensate of the form 



= - + A/J,- 



Bfi ln(/i) +C/i 2 



we indeed find 



(24) 



(25) 



The linear term in T, u (n) has cancelled, and the asymptotic behavior for ji — * gives us the constant in front of the 
/iln(/z)-term in X !/ (/i). 

For the chUE universality class, the quenched susceptibility defined above becomes quite simple: 



,ChUE 



%2^T~ = 2/ "+i (m)^-! (m) + ^ 



which has the small expansion 

^ ChUE (^) 
S 2 V 



2z/ 



i + KMf) + 7-|V 2 + o(M 4 ) 
-|(Mf) + 7V 2 + o(/i 4 ) 



1 



if u = l 

if 7/ > 2 



(26) 



(27) 



Of course, similar expressions can be derived for the partially quenched cases. 

In Fig. ||a we show raw data for u> for the same lattice couplings and lattice volumes as in Fig. [|. Again, these 
raw data beautifully collapse down on the one single scaling function u>o(a*) when rescaled according to the above 
prescription, as shown in Fig. [3|d. We emphasize again that the data for w(/i) are much more sensitive probes of the 
microscopic spectral density of the Dirac operator than the chiral condensate itself. 

For the chSE and chOE universality classes the general expressions for w(/i) are quite involved, but the cases with 
v = and 1 become relatively simple. The prediction for the chSE universality class in a sector of topological charge 
zero and one reads: 
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.ChSE/ \ 

" (: = 7r[i^o(2M) + 2Mifi(2/i)][Lo(2/x)7 1 (2/i)-L 1 (2 M )7 (2 M )]-4 M ^ 1 (2 M )7 2 (2 Al ) , (28) 



£ 2 F 

^ hSE (/i) 2 



v; , r - , , iKi(2fj) [I 3 (2n) + M I 2 (2 M )] + 8if 2 (2 M )I 2 (2 M ) 

-7T [3X 2 (2 M ) + 2/^(2^)] [L (2/i)/ 1 (2 M ) - Li (2 M )7 (2/i)] (29) 



with the small [i expansion 

, ,ChSE/ 



(p) _ ^ r-|(ln(^)+7-iV 2 + 0(M 4 ) iff = 



£ 2 F M 2 \ + C(M 4 ) ifv = l. 

The analogous prediction for the chOE case is 



(30) 



,ChOE 



0) K 



+^ [/ (m) - Mm)] - f [/iM - Mm)] + 1 , (3i) 



,ChOE 



(/*) 



E 2 ^ /i : 



= — + K (n)I 2 (n) + K^hifi) (32) 



with the small /i expansions 



^ ChOE (^) = 2u f ^ - |A* - KMf ) + 7 - + ^(M 3 ) if - = 

£ 2 F m 2 + 1 I + I(ln(f)+ 7 -i) Ai 2 + 0(M 4 ) if ^ = 1- 1 J 

We show the rescaled data for gauge group SU(2) and staggered fermions in the fundamental representation in Fig. ||c, 
and compare these rescaled data with the analytical prediction (pgj). The agreement is quite good, except for the 
very smallest lattice volume (4 4 at (3 = 1.8). Finally, in Fig. ||d we show analogous data for gauge group SU(2) and 
staggered fermions in the adjoint representation, where the analytical prediction (of the chOE universality class) is 
as given in Eq. (|3l|). Again the agreement is perfect. 

From the Leutwyler-Smilga sum rules and from the small mass expansions for the chiral condensate we can make a 
general prediction for the agreement of the condensate with the RMT predictions based on how well the microscopic 
spectral density fits the corresponding RMT predictions. We see in the unitary case that the coefficient of the linear 
order in the mass prediction for the sum rule Eq. (|l^) is dependent on all the non-zero eigenvalues (appropriately 
weighted) in the spectral sum for v > 0. For the v = case, we find the condensate will depend most strongly on the 
leading edge of the lowest eigenvalue for very small [i. This dependence is related to the appearance of the fi ln(/i) 
term in Eq. ([l0|). For the symplectic case, the coefficient in Eq. (16) is predicted for all v > (no appearance of a 



ln(/x) term at leading order). For the orthogonal case, there is no sum-rule and we can expect a strong dependence on 
the smallest non-zero eigenvalue in all topologies at small fermion mass. In general then if the microscopic spectral 
density agrees well for many oscillations with RMT we can expect reasonable agreement for the condensates with 
RMT in the ensembles and topological sectors where the sum-rules apply. In the other cases, when probing with a 
small fermion mass there can be a strong dependence on how well the smallest eigenvalue is sampled. 

At the same time, there is another competing effect that can make the non-zero topology data fit reasonably well 
with the RMT predictions. From the solution to rj(m) in Eq. (^), we see that when m — > the cutoff for the bottom 
spectrum of D^(0)I?(0) is the smallest eigenvalue which will be non-zero. Hence, r\ goes to a constant at m = 0. 
However, (tpip) and u> go down with an explicit m or m 2 factor. The RMT predictions are non-trivial because, in 
the unitary case for example, a /xln(/Lt) term is generated. At higher topology, the ln(/x) term moves to higher C(/i). 
Hence, we can see trivial agreement with RMT at higher Q when we probe the smallest eigenvalue. However, we still 
have to get the overall infinite volume scale E correct and that is non-trivial. 



III. TOPOLOGY: OVERLAP FERMIONS 



It is particularly interesting to test the analytical predictions for sectors with non-trivial topological charge. While 
staggered fermions are unsuitable for this, there now exist lattice-fermion formulations which correctly reproduce 
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FIG. 4. The microscopic spectral density for overlap fermions in (a) the fundamental representation of SU(3) (unitary), (b) 
the fundamental representation of SU(2) (orthogonal), (c) the adjoint representation of SU(2) (symplectic). Deviations from 
the analytical predictions (smooth curves) beyond a few oscillations are due to the rather small physical volumes considered. 
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FIG. 5. The distribution of the lowest eigenvalue for overlap fermions in the orthogonal ensemble. The scale is enlarged to 
show the leading edge of the distribution. 



those chiral Ward identities that are sensitive to gauge field topology. Because they share the same Ward identities 
as continuum fermions, their effective Lagrangians coincide with those of conventional chiral perturbation theory. 
In particular, in the scaling limit (1), these lattice fermions will give rise to the same Leutwyler-Smilga effective 
Lagrangians (depending on the gauge groups and color representations), and will hence fall into exactly the same 
universality classes as continuum fermions. 

The overlap Dirac operator [ ft5| derived from the overlap formalism [ pl| is a proper realization of a single flavor 
massless fermion on the lattice that separates lattice gauge fields into different topological classes based on the number 
of exact fermion zero modes. The massive overlap Dirac operator is given by 



D(m) = ~[l + m+(l-m)j 5 e(H w )} 



(34) 



with < m < 1 describing fermions with positive mass all the way from zero to infinity and where H w is the 
hermitian Wilson-Dirac operator with a negative Wilson-Dirac mass on the lattice [|l5 16 23]. Here t{x) indicates the 
sign function. 

The external fermion propagator is given by 



D~ l {m) = (1 - m)- 1 [.D -1 (m) - l] 



(35) 



The subtraction at m = is evident from the original overlap formalism [ pl| and the massless propagator anti- 
commutes with 75 |15|. With our choice of subtraction and overall normalization the propagator satisfies the relation 



m(i)| 75 ^ 1 (ra) \b) = (b\L)- 1 (m)\b) V6 satisfying 75 |6) = ±\b) 



(36) 



for all values of m in an arbitrary gauge field background |l6| ]. If chiral symmetry is broken, the right hand side of 
Eq. ( |36| ) is non-zero in the massless limit implying that the pion mass goes to zero as the square root of the fermion 
mass. Since (75-D(m)) 2 commutes with 75, its eigenvectors are chiral. In the basis where ( 7 5-D(m)) 2 is diagonal, 
7 sD(to) is block diagonal with each block being a 2 x 2 matrix (jl6). Exact zero eigenvalues of {^D{m)) 2 are paired 
with unit eigenvalues of ( 7 s.D(m)) 2 with the opposite chirality. These eigenvectors of ( 7 5D(m)) 2 are also eigenvectors 
of D(m) and therefore the topological zero modes of D(m) are chiral. We shall denote the non-zero eigenvalues of 
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(j5D(m)) 2 by A 2 with < Af < 1. In terms of these eigenvalues, the chiral condensate and chiral susceptibility in a 
fixed gauge field background are given by |]l6| 



1 v-* i ~r / \ it w 1 m 1, \v\ 1 v- * 2m(l — X 2 ) 



(1 — to 2 ) + to 2 



and 



to = -(iPi>) A - = 77 (Ml5D)- 2 (m) + TvD- 2 (m) 

m am V V 



(37) 



(38) 



Similar to Eq. (f56|), we have |3[ 

(6|( 75J D)- 2 (m)+^- 2 (m)|&) 

where 



1 — m 2 
2m 2 



(1 — m 2 )' 



(jj(rn) — b) 

(rf {m) — &') (77(771) — 6) 



H 2 (m)rj{rn) = b with H (m) = 7 5 D(m), 756 = ±6 



We note that 



with 



H 2 (m) = D\m)D{m) = D{m)D\m) = (l - m 2 ) 



#o(0) 



1 — rrr 



^ o 2 (0)Xd 



^ + ^(75±l))e(^) 



X± with 75X± = ±X± 



(39) 



(40) 



(41) 



(42) 



Eq. ([ll]) implies that we can solve the set of equations H 2 (m)r)(m) = b for several masses, to, simultaneously (for the 
same right hand b) using the multiple Krylov space solver described in Ref. |20) . In our tests we used fermion masses 
from to 0.999. However, for comparisons to RMT we only consider the fermion mass range from to = 10~ 4 up to 
0.999. 

The first term on the right hand side of Eq. (|3^) is due to the presence oi\v\ exact zero modes in a fixed gauge field 
background. By working in the chiral sector where {^^D{m)) 2 has no zero modes, it is possible to compute the second 
term in Eq. (|7|) and investigate the onset of chiral symmetry breaking on the lattice Note that the relation 

( |37| ) is an exact identity at any lattice spacing. It is of the same form as Eq. (^|), up to terms vanishing with the 
ultraviolet cut-off. The bare fermion mass enters the overlap Dirac operator in a non-trivial way and is proportional 
to the mass parameter to in ( |34"| ) only for small to, i.e. only up to terms of relative 0{a 2 ). The proportionality factor, 
Z m , depends in particular on the mass in the Wilson-Dirac operator used |l6f| . Since Z m is the inverse of the wave 
function renormalization constant |l6| ], the rescaled mass parameter fi = mT,V is independent of these Z-factors and 
agrees with the continuum definition up to terms vanishing with the ultraviolet cut-off. 

The infinite- volume chiral condensate S differs significantly, at the same /3- values, from that of staggered fermions. 
However, in the cases we shall present here this one single parameter E has already been extracted to high precision 
from the distribution of the smallest Dirac operator eigenvalue jT3 ], The analytical predictions for are thus 

parameter-free also in these cases. 

On all the gauge configurations used for the measurement of the condensate a few low lying eigenvalues have 
previously been determined We thus know the number and chirality of all zero modes, and hence the topological 
charge. As mentioned already, when zero modes are present, we perform the stochastic estimate in the sector with 
opposite chirality. In topologically trivial gauge fields, we perform the stochastic estimate in the positive chirality 
sector. 

We are now ready to test some of the predictions for overlap fermions in the finite-volume regime. The first 
observation is that the universality classes of continuum fermions coincide with those of overlap fermions. We begin 
by comparing the microscopic spectral density p s (C) with the predictions of RMT. In Fig. |^ are the results for different 
topological sectors for all the ensembles. The curves are the predictions using the infinite volume S previously 
determined (T^|. We see good agreement for the first oscillation (essentially the lowest eigenvalue contribution) for 
all topological sectors and the best agreement in the symplectic case in Fig. 0c. However, the data rapidly deviates 
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FIG. 6. The rescaled condensate for overlap fermions in (a) the fundamental representation of SU(3) as function of [s, = m'EV, 
(b) the fundamental representation of SU(2), (c) the adjoint representation of SU(2). 
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up from the curve for higher eigenvalues. This is a finite- volume effect since as the volume increases the scale of the 
eigenvalues decreases into the region where RMT applies. The infinite volume chiral condensate S sets the scale for the 
eigenvalues, and since the condensate is larger for staggered fermions compared to overlap fermions at corresponding 
parameters, the eigenvalues for staggered fermions occur closer to zero where there is a corresponding better agreement 
for more oscillations with the RMT predictions. 

Nevertheless, we can expect to find the best agreement of lo v with the RMT predictions where the fermion mass is 
probing the scale of eigenvalues in p s {C) that are in best agreement. For the chiral condensate (ipip)i' the agreement 
with RMT will depend on the overall p s {C) which we remarked is best for the symplectic case. Since the spectral 
density does not match the RMT prediction for many oscillations in the unitary case as seen in Fig. ||a, we expect 
higher topologies to not be well described by RMT. 

In addition, we can expect problems for comparisons of observables for small fermions masses when the p s (z) has not 
been adequately sampled, and this problem is most pronounced in the orthogonal case. Compared to the staggered 
fermion example in Fig. |^, we do not have enough statistics and large enough volumes to adequately sample the 
leading edge of the eigenvalue distributions shown in detail in Fig. [|. 

We continue with tests of the predictions for E„(/i) in the chUE case, using quenched overlap fermions and gauge 
group SU(3). Shown in Fig. ||a are some data for gauge field sectors with v = . . .3. We stress that for the sectors 
of non-vanishing v we have subtracted the somewhat trivial v / p term, which otherwise would completely dominate 
the plot. What is shown is thus not the chiral condensate per se, but rather E„(/i)/£ — |f The agreement in the 
v = sector is good, but while the data for the v = 1, 2, and 3 qualitatively display the right behavior, they are 
nevertheless somewhat off the analytical predictions. 

In Fig. ^a we show the subtracted u v /T?V for the unitary case. As mentioned before, this observable is a more 
sensitive probe of the eigenvalue distribution at a fixed fermion mass and should yield better agreement when the 
corresponding p s (C) is in good agreement with the RMT predictions. We see good agreement for z = m^V > 0.05 
in all topology sectors, but below this value there are deviations related to the lack of small eigenvalues in p a and 
the small volume. For the orthogonal ensemble we find even worse agreement for £„(/i)/£ — \v\j p shown in Fig. ||b, 
but find a small mass region (0.4 < z < 4) in ui v jY?V shown in Fig. [7|b where we simultaneously have an adequately 
sampled distribution of very small eigenvalues and good agreement of p s (z) with RMT. 

We next turn to gauge group SU(2) and overlap fermions in the adjoint representation where we find the expected 
good agreement with analytical predictions for the v = and 1 sectors (we found almost no configurations in the 
sectors of higher topological charge in this case). These graphs are shown in Fig. ||c and Fig. 0c. 



IV. CONCLUSIONS 



We have performed a systematic series of Monte Carlo tests of the analytical predictions for the chiral condensate 
and related chiral susceptibilities in the finite- volume scaling region of Eq. (1). In four dimensions there are three 
universality classes with which to compare, conveniently classified in Random Matrix Theory terminology by means 
of chiral versions of the three classical matrix ensembles, i.e., chSE, chUE and chOE. Once the infinite- volume chiral 
condensate £ is known, there are parameter-free finite-volume scaling functions with which to compare data. As we 
have shown, results for all three universality classes with topological charge v = are nicely reproduced by staggered 
fermions. To test the analytical predictions for gauge field sectors of non-trivial topological winding numbers, we 
have also considered overlap fermions, which possess exact zero modes in topologically non-trivial gauge fields. Here 
there is qualitatively good agreement, with even excellent agreement in the case of the chSE universality class. The 
deviations from the RMT predictions observed in the other cases are understood to be due to either finite volume 
effects or the relatively modest statistics we were able to obtain in the simulations with overlap fermions. 

The results presented here clearly show the power of the finite-size analysis that has come out of the study of 
finite-volume effective partition functions and Random Matrix Theory. In contrast to conventional finite-size scaling 
analysis near critical points, we are here in the unusual situation of knowing not only the right scaling variables, 
but also parameter-free exact analytical predictions for the scaling quantities. In this particular corner of those 
non-Abelian or Abelian gauge theories that support spontaneous breaking of chiral symmetry the exact analytical 
predictions have thus very clearly been confirmed by direct numerical studies. 
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